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Abstract. In the supersymmetric quantum mechanics formalism, the shape 
invariance condition provides a sufficient constraint to make a quantum mechanical 
problem solvable; i.e., we can determine its eigenvalues and eigenfunctions 
algebraically. Since shape invariance relates superpotentials and their derivatives at 
two different values of the parameter a, it is a non-local condition in the coordinate- 
parameter (x, a) space. We transform the shape invariance condition for additive shape 
invariant superpotentials into two local partial differential equations. One of these 
equations is equivalent to the one-dimensional Euler equation expressing momentum 
conservation for inviscid fluid flow. The second equation provides the constraint that 
helps us determine unique solutions. We solve these equations to generate the set 
of all known ?i- independent shape invariant superpotentials and show that there are 
no others. We then develop an algorithm for generating additive shape invariant 
superpotentials including those that depend on h explicitly, and derive a new h- 
dependent superpotential by expanding a Scarf superpotential. 
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1. Introduction 

1.1. Background 

Supersymmetric quantum mechanics (SUSYQM) provides an elegant and useful 
prescription for obtaining closed expressions for the energy eigenvalues and 
eigenfunctions of many one dimensional problems and three dimensional problems with 
spherical symmetry [H Ej [3l H] . 

In the SUSYQM formalism, first order differential operators A~ and A'^ are 
generalizations of the raising and lowering operators employed by Dirac for treating 
the harmonic oscillator [5J. These "ladder" operators y4^(x, oq) = =F^^ + Wi^x^ao) use 
a superpotential W to generate partner Hamiltonians H_ = A'^A'^ and = A^A'^. If 
these partner potentials are shape invariant, then the eigenspectra for both Hamiltonians 
can be derived algebraically without any prior information for either Hamiltonian. 

Until recently, all known shape-invariant potentials could be generated from 
superpotentials with no explict dependence on ^ [31 H] . We classify such superpotentials 
as "conventional." However, a new class of shape-invariant potentials was discovered by 
Quesne [6] and expanded elsewhere [3 E\- These potentials arise from an "expanded" 
set of superpotentials that contain explicit ^-dependence. 

In a recent publication [9], we showed that the shape- invariance condition can be 
transformed into two local partial differential equations. Solutions to these equations 
generate the set of all known conventional shape invariant superpotentials and allow no 
others in this category. In addition, these equations provide an algorithm for generating 
^-dependent potentials. In this manuscript, we elaborate this method, proving the 
completeness of the set of "conventional" superpotentials and generating a previously 
unknown "expanded" superpotential. 

1.2. Supersymmetric Quantum Mechanics 

A quantum mechanical system in one spatial dimension x described by a potential V{x) 
can alternately be described by its ground state wavefunction ^o- For simplicity of 
notation, we will use units where 2m = 1 throughout this manuscript. We may also set 
the ground state energy to zero without changing the physics. With these choices, the 
Schrodinger equation for the ground state wavefunction is 



where prime denotes differentiation with respect to x. Thus, it follows that the potential 
can be written as 



The SUSYQM formalism makes use of first order differential "ladder" operators 
A~ and A+, 



n% + v{x)^o = , 



(1) 




(2) 




(3) 
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These ladder operators include a superpotential Oq), which is a real function of 
coordinate x and parameter (or a set of parameters). 

Operators generate two supersymmetric partner Hamiltonians: = A^ A~ 
and if+ = A~ A^ . Hamiltonian if+ is called the superpartner of H_, and corresponding 
potentials V- and V+ are given by V± = W'^{x^ oq) ± • 

Let us denote the eigenf unctions of H± that correspond to eigenvalues E^, by ^Z'^^-'. 
For positive integer n, 

{A-4;)) = A-A^ 

= A- 

Hence, with the exception of the the ground state which obeys A'lpl ■* = 0, all excited 
states -ip!^^ of H_ have one-to-one correspondence with eigenstates of with exactly 

), as illustrated in 
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Figure 1. Schematic illustrating the isospectrality of and showing the ground 
state and first two excited states of i7_ (corresponding to potential V-{x) sketched 
with a solid bold curve) in the left column, and the ground state and first excited state 
of i/+ (corresponding to potential V+(a;), sketched as a sohd bold curve) in the right 
column. Sample energy levels are shown as horizontal lines, and a sketch of a sample 
wavefunction is overlaid with a dashed line for each energy level. In general, V+ and 
have different shapes, as do various V'^ and ^~ . However, the energy = Eq , 
E2 = E^ , and E^_-^ = E^ . The operator A~ acting on an eigenstate of will yield 
an eigenstate of , and acting on an eigenstate of will yield an eigenstate of 
H+. 

Conversely, oc A'^ipl^-^. Thus, if the eigenvalues and the eigenfunctions of H_ 
were known, one would automatically obtain the eigenvalues and the eigenfunctions of 
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(and vice versa), which is in general a different Hamiltonian. However, unless we 
know one set of eigenstates a priori, this analysis is simply a mathematical curiosity. 

1.3. Shape Invariance 

Despite this limitation, for special cases in which these partner potentials obey the 
"shape invariance" condition [101 E] • 

V+{x, flo) + g{ao) = Vl(x, ai) + g{ai) , (5) 

where parameter ai is a function of oq, i.e., ai = /(ao); the spectra for both Hamiltonians 
can be derived algebraically without any prior information for either Hamiltonian. This 
is due to the existence of an underlying potential algebra [121 [131 E] • 

Let us consider the special case where V^_(x, oq) is a shape invariant potential. 
In this case, potentials V- and V+ have the same x-dependence, and the corresponding 
Hamiltonians H+{x, Oq) and H_{x, ai) differ by g{ai) —g{ao). Thus, their eigenfunctions 
are the same, and corresponding eigenvalues differ by g{ai) — g{ao). In particular, 
they have a common ground state wavefunction, given by iljQ^\x,ao) = iljQ~\x,ai) ~ 
exp (^— J^^W{x,ai)dx^, and the ground state energy of if+(x,ao) is g{ai) — g{ao), 
because the ground state energy of H^{x, ai) is zero. Note that the parameter shift oq — >■ 
ai has an effect similar to that of a ladder operator: iIj[~\x, Qq) ~ v4+(x, ag) tljQ~\x, ai). 
Also note that the ladder operators A~ and A'^, like H^, are dependent on parameters 

The first excited state of H^{x,ao) is given to within normalization by 
A'^{x,aQ)iljQ~\x,ai) and the corresponding eigenvalue is g{ai) — g{aQ). By iterating 
this procedure, the (n + l)-th excited state is given by 

V'i+\(x,ao) ~ A+(ao) A+(ai) ■ ■ ■ A+(a„) 4~^(a;,a„) , (6) 

and corresponding eigenvalues are given by 

Eq = 0; and E^") = g{an) - g{ao) for n > 0. (7) 

(To avoid notational complexity, we have suppressed the a;-dependence of operators 
A{x, ao) and A~^{x, Qq).) Thus, for a shape invariant potential, one can obtain the entire 
spectrum of if_ itself by the algebraic methods of SUSYQM (and of course the same is 
true for H^). 

In this manuscript, we develop a method for finding shape- invariant superpotentials 
from a system of partial differential equations, and use this method to discover a new 
shape-invariant superpotential. In Section II, we will show that for /^-independent 
("conventional") superpotentials, the shape invariance condition can be converted into 
an infinite sequence of partial differential equations. In Section HI, we solve these partial 
differential equations and systematically generate the complete set of conventional shape 
invariant potentials. In Section IV, we show a method for using these equations to find 
/i-dependent ("extended") superpotentials. In Section V we use this method to find a 
new shape-invariant superpotential, and we present our conclusions in Section VI. 
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2. Expressing the Shape Invariance Condition With Partial Differential 
Equations for Conventional Superpotentials 

As we have just seen, any partner potentials obeying ([5]) can be solved algebraically. 
Thus, discovering that a potential is shape-invariant yields much useful information. 
Here we develop a method for finding shape-invariant potentials. 

For shape invariant systems, the energy eigenvalues of -ff_(x, Oq) are given by 
E^~\aQ) = g{an) — g{ao), where a„ = f"'{ao) indicates / applied n-times to oq [3111]. 
To avoid level crossing, g{a) must be a monotonically increasing function of a; i.e., 
1^ > 0. In the case of "additive" or "translational" shape invariance, the parameters 
differ by an additive constant; i.e., aj+i = ai + h. While there are other forms of shape 
invariance such as multiplicative [15] and cyclic [16], most of the known exactly solvable 
superpotentials exhibit additive shape invariance [IHIITT]. We thus restrict our analysis 
to additivel shape invariance. 

Several groups found these potentials by imposing various ansatzes [TBI [HI [I9]. In 
this manuscript, we derive these potentials ab initio as the solutions to a set of partial 
differential equations that must be satisfied for all additive shape-invariant potentials. 
Writing (jS]) in terms of the superpotential, we obtain 

„„) + n + <,(a„) = W\., a,)-n '^^^ H- 9(a,). (8) 

Note that ([8]) is a difference-differential equation; that is, it relates the square of the 
superpotential W and its spatial derivative computed at two different parameter values: 
(x, Co = a) and (x, ai = a + h). This equation has also been studied for dynamical 
systems, where it is known as the infinite-dimensional dressing chain [20| |2T] . Note that 
we have not needed to specify the value of h\ therefore, this equation must hold for 
any value. Using this feature, we will transform (|8]) into a non-linear partial differential 
equation that is local in nature; i.e., all terms can be computed at the same point (x, a). 
This has the obvious advantage of mathematical familiarity (at least to physicists). In 
addition, it provides a systematic method for finding additive shape invariant potentials. 

Important correspondences have been shown to exist between quantum mechanics 
and fiuid mechanics [22]. SUSYQM is well known to have a deep connection with 
the KdV equation [231 EH ESI [271 lEE], a nonlinear equation that describes waves in 
shallow water. In this section we show that every additive shape invariant superpotential 
that does not depend explicitly on h corresponds to a solution of the Euler equation 
expressing momentum conservation for inviscid fiuid fiow in one spatial dimension. We 
use this correspondence to develop a systematic method which 1) yields all such known 
/i-independent solvable potentials for SUSYQM and 2) shows that no others exist. We 
will consider ^-dependent potentials in Section III. 

Because of additive shape invariance, the dependence of on a and % is through 
the linear combination a + therefore, the derivatives of W with respect to a and h 
are related by: = dw{x^a+h) ^ gj^^^g q must hold for an arbitrary value of /l, if 

we assume that W does not depend explicitly on h, we can expand in powers of % and 
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the coefficient of each power must separately vanish. Expanding the right hand side in 
powers of h yields 

dW dW 1 dgia) 

W \ = f9l 

dx ^ 2 da ' 

dW 1 dg(a)\ , , 

Thus, all conventional additive shape invariant superpotentials are solutions of the above 
set of non-linear partial differential equations. Although this represents an infinite set, 
note that if equations at and 0{^^ are satisfied, all others automatically follow. 
Replacing W by — u, x by t, and a by a; in ([9]), the equation then becomes: 

This equation is equivalent to the equation for inviscid fluid flow in the absence of an 

external force on the bulk of the fluid: 

(9u(x, t) , , „ , , Vp(x, t) 

\: ' + « X, t) ■ Vu X, t) = 13 

ot p (x, t) 

in one spatial dimension with the correspondence = where u is the fluid velocity 
at location x and time t, p is the pressure, and p is the local fluid density. Equation f[T^ 
is one of the fundamental laws of fluid dynamics, and was first obtained by Euler in 
1755 [28j. Thus, all conventional shape invariant superpotentials form a set of solutions 
to the one-dimensional Euler equation. 

It should be noted that ( fT3i) is not solvable unless additional constraints are 
applied. In fluid dynamics this equation is generally supplemented by the continuity 
equation expressing conservation of mass, along with an equation of state and/or the 
energy equation and boundary conditions. These additional constraints do not apply 
in SUSYQM. Instead, (|TT1) supplies the additional constraint which must be fulffiled to 
satisfy shape-invariance. Thus, the set of solutions of (jH]) that also satisfy the constraint 
of (fTTj) will define the complete set of conventional shape invariant superpotentials. 

3. Generating the complete set of conventional superpotentials 
3.1. Solutions of Special Cases 

In this section, we show that the set of all possible conventional superpotentials are 
determined by six special cases. To find this set of solutions, we note that ffTTl) is 
satisfied for all n > 3 as long as 

The general solution to Eq. ( IT^ is 

W{x,a) = a- Xi{x) + X2{x) +u{a) . (15) 



Method for Generating Additive Shape Invariant Potentials from an Euler Equation 7 



Substituting this into ^ yields 

^ du\ ( dXi dX2\ I do 

V ' Y da J \ dx dx J 2 da 

W ^ ^ ' ^ V ' 

aw dw 

da dx 

To systematically find all possible solutions, we define m ^ + | ^ = ~H{a), and then 
collect and label terms based on their dependence on Xi and X2 and their derivatives: 

Term#l ^ ^ ' Term#3 ^ « ' t^~^c; ^ ^ ' 

Term#2 ^ Term#4 Term#5 Term#6 

We begin by considering special cases of (fT7|) where one or more of the terms Xi{x), 
X2{x), or u can be considered zero. After considering these cases, we will show that all 
solutions to this equation can be reduced to one of these cases. In the nomenclature 
that follows, lower case Greek letters denote a- and x-independent constants. 

3.1.1. Case 1: X2 and u are not constants, Xi is constant. In this case, let Xi = fi. 
Then our general form for W becomes W = fia+u{a)+X2{x). If we define u = u{a)+fia, 
we get W = U + X2. So this case is equivalent to Xi = 0. Then terms 1, 3, 4, and 6 each 
become zero, and (fTTll becomes — + = H(a). Since X2 must be independent 
of a but cannot be constant,^- 7^ 0. Thus, this is possible only if ^ and H{a) are 
independent of a. This yields u = aa + 13, H{a) = 9. Therefore, — + = 9, 
where a 7^ since u is not constant. The solution to this differential equation is 
-^2(3;) = ~ + i]e"^. Therefore, W = aa + (3 + ^ + r]e°'^. Defining a = —1, this yields 
W = A — Be^^, where A = P — a — 9. This is the Morse superpotential. 

3.1.2. Case 2: Xi and u are not constants, X2 is constant. Following a similar 
procedure this case is equivalent to X2 = 0. So aXf — + (u + Xi = H{a). 
Since Xi cannot depend on a but must contain a;-dependence, the only way for H{a) to 
be independent of a: is if u + a^ = aa, where a could be any constant. Thus, u = ^ + f , 
where a and (3 could be any two constants, although they cannot both be zero. 

Thus, aX| — 0,^^ + aaXi = H(a). This is only possible if H{a) = 7a for some 
constant 7. So Xf — + aXi = 7. This differential equation gives different solutions 
depending on the values of a and 7. 

If a 7^ 0, then 7 7^ yields the Rosen-Morse I superpotential, and 7 = yields the 
Eckart superpotential. 

If a = 0, then 7 = yields Coulomb, 7 > yields Rosen-Morse II, and 7 < yields 
Rosen-Morse I. 

Thus, the Rosen-Morse I, Rosen-Morse II, Eckart, and Coulomb superpotentials 
are all solutions to Case 2, for different values of a and 7. 
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3.1.3. Case 3: Xi and X2 are not constants, u = + v (this includes the case where 
u is constant.) In this case we can define Xi = Xi + ^ and X2 = X2 + z/, which is 
equivalent to m = 0. For m = 0, X1X2 - ^ + a (x^ - ^) = H{a). 

Since Xi and X2 are both independent of a, the coefficients of each power of a must 
cancel separately. So we are left with two coupled differential equations: Xf — = a, 
X1X2 — = /3. Again, the solution varies depending on the values of a and /3, which 
could be any constants. 

If a = 0, the solution is the 3-D harmonic oscillator superpotential. If a < 0, the 
solution is the Scarf I superpotential, and if a > 0, the equation is solved by either the 
Scarf II or Generalized Poschl- Teller superpotential. 

Therefore, the Scarf I, Scarf II, 3-D oscillator, and Generalized Poschl- Teller 
superpotentials are all solutions of Case 3. 

3.1.4. Case 4- X2 is not constant, Xi and u are constant. If Xi 7^ 0, this is equivalent 
to Xi = and u = fia + u. In such a case, this is equivalent to Case 1, and the solution 
is the Morse superpotential. However, if Xi = 0, this case is equivalent to Xi = u = 0, 
in which case — = H{a). Since X2 is independent of a, H{a) must be a constant. 
This generates the one- dimensional harmonic oscillator. 

3.1.5. Case 5: Xi is not constant, X2 and u are constant. In this case, aXi + aXf — 
a^- = H{a) for some constant a. Since Xi is independent of a but must depend on x, 
a = and H{a) = /3a for some constant (3. Thus, Xf — = (3. This yields special 
cases of Scarf I and Scarf II, and the centrifugal term of the Coulomb and 3-D oscillator, 
depending on whether (3 is positive, negative, or zero. 

3.1.6. Case 6: Xi is constant, X2 is constant. In this case, the superpotential has no 
x-dependence, regardless of the value of u. This is a trivial solution corresponding to a 
fiat potential, and we disregard it. 

Thus far, we have considered all of the special cases that are equivalent to m = 0, 
Xi = 0, or X2 = 0. These special cases generate all known conventional additive 
shape- invariant superpotentials [3[ [T71 H] . We have listed them all in Table [H 

3.2. Proof That the List of Conventional Superpotentials is Complete 

Now that we have these special cases, we can systematically obtain all possible solutions. 
H{a) is independent of x. Therefore, when any solution is substituted into (fT7|l . it will 
yield an x-independent sum of terms 1-6. There are many ways in which these terms 
could add up to a term independent of x. As a first step, we begin with the simplest 
possibility, in which each term is individually independent of x. 

Under this assumption, term 3 states that Xi must be a constant, independent of 
X. In addition, term 1 dictates X1X2 must be constant as well. These two statements 
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Name 


superpotential 


Special 

Cases 


Harmonic Oscillator 


2 


Xi = M = 


Coulomb 


2{£+l) r 


X2 = 


3-D oscillator 


2 r 


■u = 


Morse 




Xi = 


Rosen-Morse I 


-A cot X - f 


X2 = 


Rosen-Morse II 


A tanh x f 


X2 = 


Eckart 


-Acothx -f f 


X2 = 


Scarf I 


Atanx — i?secx 


M = 


Scarf II 


A tanh x + i?sech x 


M = 


Gen. Poschl- Teller 


y4 coth X — i?cosech x 


M = 



Table 1. The complete family of additive shape-invariant superpotentials. 



can only be true if X2 and Xi are constant separately; this reduces to the trivial solution 
of case 6. 

Therefore, assuming that each term is separately independent of x yields only the 
trivial solution. However, there is also the possibility that some of the terms depend 
on X, but when added to other terms, the x-dependence cancels to yield a sum that is 
independent of x. If a group of ra-terms taken together produce an x-independent sum, 
and if no smaller subset of these terms add up to a sum independent of x, we call this 
group of n-terms "irreducibly independent of x" . 

As an example, we check whether there are any solutions for the six-term irreducible 
set in which the sum of all six terms in ( ITTl) is independent of x, but in which the sum 
of any subset of terms would depend on x. 

To check this possibility, we note that the first two terms are independent of a, while 
terms 3 and 4 are linear in a. We do not know a priori the functional form of u. Since 
term 5 contains it could include terms independent of a, terms linear in a, and/or 

da ' 



other forms of a-dependence. Similarly, we do not know the functional form of m + a-'^ 



Therefore, we define functions /i(a) and /2(a) such that u + = ai + j3ia + /i(a) and 
^ = Ci2 + 1^20- + /2(a), where we do not know the functional form of /i(a) and /2(a), 
except that they contain no constant terms or terms linear in a. 

With this definition, term 6 becomes aiXi + a/3iXi + /i(a)Xi and term 5 becomes 
02X2 -|- 0/32X2 -|- /2(a)X2 Since /i(a) and /2(ci) contain no constant terms or terms 
proportional to a, the x dependence of the term /i(a)Xi cannot be cancelled by 
terms 1-4, and neither can the x-dependence of /2(a)X2. Therefore, we conclude that 
/iHXi + /2(a)X2 = /i. 

This leaves only two possibilities. First we consider the possibility that = 
/2(a) = 0. In this case, 4^ = 0^2 + 1^20-^ so u = + + 7. Plugging this solution into 



term 6 then becomes (^a^ + 2a;2 + 7) Xi. Since Xi is independent of a and none 
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of the other terms are proportional to a^, must be independent of x. So either 

Xi is a constant, or = 0, in which case u is hnear in a; either possibihty represents a 
special case. 

On the other hand, if 7^ 0, then either Xi is a constant (which is a special case), 
or the x-dependence of /i(a)Xi must be cancelled by the x-dependence of /2(a)X2. Since 
Xi and X2 cannot depend on a, this is only possible if X2 + /iXi = v for some constants 
/i and V. However, if this is the case, then W = aXi + X2 + u = aXi — fiXi + u + u. 
Since the zero of a can be shifted by defining a = a — ^, and we can define u = u + u, 
this is equivalent to the case W = dXi+u. So this particular case can be reduced to the 
special case of X2 = and does not yield any new solutions. Thus, no new solutions can 
be found by assuming the full six-term set is irreducibly independent of x; any solution 
found under this assumption can be reduced to one of the special cases. 

Since each term cannot be separately independent of x and the full six-term 
set cannot be irreducibly independent of x, the only possibility that could produce 
solutions not covered under one of the special cases would be if two or more irreducibly 
independent sets combine to produce a solution. If, for example, term 5 depends on x 
and term 6 depends on x, but the sum of these two terms is x- independent, then we 
consider the set of terms {5, 6} to be a two-term set that is irreducibly independent 
of X. To continue this example, if the four sets {1}, {2,3}, {4}, and {5,6} were each 
irreducibly independent of x, then the entire left side of (|T71) would be independent of 
X, and this could possibly produce a solution not covered under the special cases. We 
now proceed to examine all possible combinations and show that no new solutions are, 
in fact, produced. 

To examine all possibilities, we begin by examining all possible one-, two-, and 
three- term sets, and see what restrictions are imposed on solutions in each case. For 
example, if term 1 is individually independent of x, then X1X2 = a for some constant 
a, which implies that Xi = a/X2. Table [2] shows the consequences required for each 
possible one-term set to be independent of x. 

We now proceed to examine two-term sets. Let us take the example mentioned 
above in which {5, 6} is a two-term set that is irreducibly independent of x (that is, 
term 5 depends on x and term 6 depends on x, but the sum of these two terms is 
x- independent). Let us consider this example further to see if it leads to any new 
solutions. In this example, ^-^2 + (u + Xi must be independent of x. Since Xi 
and X2 must each depend on x and ^ must depend on a (or this would reduce to a 
special case), the only way this is possible is if the x-dependence of ^-^2 is cancelled 
by the x-dependence of (^u + Xi. For this to be true, the a-dependence of ^ must 
differ from the a-dependence of u + by only a multiplicative constant. Thus, we 
conclude that X2+aXi = /3 for some constants a and /3. As we did when considering the 
full six-term set, we can absorb X2 into Xi by shifting the zero of a. So this particular 
case can be reduced to the special case of X2 = and does not yield any new solutions. 

However, not all two-term sets can be reduced in this way. As an example, we 
consider the set {1,2}. For this set to be irreducibly independent of x, X1X2 and 
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GoniDatible with 

* — v^xxx caj u X x\^ V V X u XX 


set 


x-independence 


a new solution? 


|1| 

L J 


Xi = a 1X2 


Ifa^O 


{2} 


X2 = ax + f3 


If a ^ 


{3} 


Xi = a 


No 


{4} 


Xi = ax + f3 


If a ^ 


{5} 


X2 = a 


No 


{6} 


u = a/a 


If a ^ 



Table 2. Assuming that a single term is independent of x results in restrictions 
on the possible solutions. The first column lists each possible single term set. The 
middle column of each row shows what consequences are necessary for that set to be 
independent of x. Throughout this section, lower-case greek letters indicate constants 
that are independent of both a and x. Finally, if this requirement is compatible with 
solutions not included as one of the special cases, restrictions on such a solution are 
shown in the third column of that row. The third column states "No" if the requirement 
is equivalent to a special case and thus no new solutions are allowed. 

must each depend on x. However, X1X2 — must be independent of x. Since Xi and 
X2 do not depend on a, the only way for this to be possible is if X1X2 — = a for 
some constant a. By itself, this requirement could allow many solutions. Therefore, we 
will have to check whether this set can couple with the remaining terms 3, 4, 5, and 6 in 
such a way as to produce new solutions that are not compatible with the special cases. 

In Table [21 we display the consequences required for each possible two-term set 
to be irreducibly independent of x. From this table, it is clear that there are several 
two-term sets that can be irreducibly independent of x. However, in order to satisfy 
(fTTj) . these sets must be compatible with solutions that allow the remaining four terms 
to be independent of x. We now check whether there are any combinations of one- and 
two- term sets that can satisy (|T71) . Once we have completed this, we will consider cases 
involving three-, four-, and five- term irreducible sets. 

Combining the results from Table [2] and Table [3] allow us to immediately eliminate 
many possibilities. For instance, in the example listed above, the set {1, 2} is irreducibly 
independent if X1X2 — = a for some constant a. However, in order to satisfy (IT7|) . 
the remaining terms must combine in such a way that the combination of terms 3, 4, 5, 
and 6 are independent of x as well. We first note from Table [2] that term 3 cannot be 
independent of x by itself, and neither can term 5. From Table |3l we note that {3, 6} 
cannot be irreducibly independent of x, and the set {5, 6} is equivalent to a special case. 

Therefore, there are only two possible combinations of one- and two- term sets 
involving the irreducible set {1, 2} that could satisfy f|T7j) and lead to a new solution. The 
first possibility is that set {1,2} is irreducibly independent of x, set {3, 5} is irreducibly 
independent of x, and set {4, 6} is irreducibly independent of x. However, from Table 
m {3, 5} requires u = ^ + 7 with « 7^ 0. On the other hand, {4, 6} requires u = f + 7, 
so these two terms are incompatible. The only other possibility is that set {1,2} is 
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Two-term 


Requirements for 


Compatible with 


set 


irreducible x-independence 


a new solution? 


|1,2| 


1 ^ dx " 


For any a 




Contradiction: must reduce 


No 




Contradiction: must reduce 


No 


11 51 


11 = nd + z/ 


No 


{1,6} 


« = f + 7; 

yA 1 TT 


If a ^ 0; /3 ^ 




Contradiction' must reduce 


No 


{2,4} 


Contradiction: must reduce 


No 


12 51 


11 = ///I -J- /y 


No 


12 61 


M = - + ^• 


If a 7^ 


1'-', 


- ^ = a 

"^^1 dx " 


For any 


{3,5} 


M = ^ + 7; 

X? + aXo = B 


If a ^ 


13 61 


Cyontradiction' must reduce 

* — v^xx vj X cajV-XX v_> u xv_yxx ■ xxx i^t- u x v. a i^i. 


No 


{4,5} 




If a ^ 


{4,6} 


^ = f + 7; 

Xi = /3e°^ + 


If a ^ 0; /3 ^ 


{5,6} 


X2 = aXi + /3 


No 



Table 3. Assuming that a two-term set is irreducibly independent of x results in 
restrictions on the possible solutions. The first column lists each possible two-term 
set. The middle column of each row shows what consequences are necessary for that 
set to be irreducibly independent of a;. The middle column states "Contradiction: 
must reduce" if the assumption that the two-term set is irreducibly independent of 
X leads to the contradictory conclusion that the set must be reducible. Finally, if 
this requirement is compatible with solutions not included as one of the special cases, 
restrictions on such a solution are shown in the third column of that row. The third 
column states "No" if no new solutions are allowed, either because the assumption 
of irreducibility leads to a contradiction or because the requirement on solutions is 
equivalent to a special case. 

irreducibly independent of x, set {3, 5} is irreducibly independent of x, and sets {4} 
and {6} are each separately independent of x. However, once again the requirement 
u = a/a for term 6 to be x- independent is incompatible with the requirement u = ^-1-7 
with a ^ for {3, 5}. Therefore, there can be no possible new solutions resulting from 
combinations of one- and two- term sets involving the irreducible set {1,2}. 

Comparing Table [2] and Table [3l every possible combination of one- and two- term 
sets leads to one of the following three results: a) one of the terms is directly equivalent 
to a special case (such as any combination involving the single-term set {3}); b) one of 
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the terms leads to a contradiction in which a term assumed to be irreducible must be 
reducible (such as any combination involving the two-term set {1, 3}); or c) two elements 
in the combination require different functional forms of m, leading to the impossibility 
of a common solution (such as any combination involving the sets {3,5} and {4,6}). 
Therefore, we conclude that no new solutions can be found from combinations of one- 
and two- term sets. 

However, there is still the possibility that there are solutions provided by 
combinations involving three-, four-, or five-term irreducible sets. We now test these 
possibilities, begining with three-term sets. In Table HJ we display the consequences 
required for each possible three-term set to be irreducibly independent of x. 

Following the same procedure as for two-term sets, we see that once again all 
combinations produce one of the following results: a) one of the terms is directly 
equivalent to a special case; b) one of the terms leads to a contradiction in which a 
term assumed to be irreducible must be reducible; or c) two elements in the combination 
require different functional forms of m, leading to the impossibility of a common solution. 

We now turn to the case of four-term irreducible sets. In this case, rather than 
calculating the restrictions on all possible four-term sets, we use the fact that each four- 
term set must combine with either a two-term set or a pair of single-term sets. We can 
use this fact to eliminate many possibilities. For instance, in order for set {2,4,5,6} 
to yield new solutions, it must combine either with the two-term set {1, 3} or with the 
pair of single term sets {1} and {3}. Since neither possibility can yield new solutions, 
we can eliminate the four-term set set {2,4,5,6} from consideration. 

Of the remaining four-term sets, many more can be considered by assuming 
the requirements from the complementary one- or two- term sets and plugging these 
solutions in to the remaining equation. For instance, since term 5 must depend on x 
(c/ Table [2]), the set {1,2,3,6} can only be irreducibly independent of x if {4,5} is 
irreducibly independent of x as well. However, from Table [3l this would require u = 
^ +7. Plugging this solution into (117]) requires that X1X2 - ^ +7X1 +aX2 + 
be independent of x in order for {1, 2, 3, 6} to be irreducibly independent of x. However, 
as term 3 is the only term linear in a, the only solution is for Xi to be constant, which 
leads to special case 1 and contradicts our assumptions. Therefore, {1,2,3,6} cannot 
be irreducibly independent of x. 

Finally, there are some four-term sets that cannot be eliminated in this manner. 
For example, {1, 2, 5, 6} could be irreducibly independent of x if a solution can be found 
that is compatible with the requirement that Xf — is constant which results from the 
two-term set {3,4}. To check this possibility, we assume that {1,2,5,6} is irreducibly 
independent of x. In this case, X1X2 — + %X2 + [u + Xi must be independent 
of X. We do not know a priori the functional form of u{a) or but we do know that 
^ must depend on a or this would reduce to a special case. Since terms 1 and 2 are 
independent of a, the only way for this to be the case is if the x-dependence of ^X2 is 
cancelled by the x-dependence of {u + Xi. For this to be true, the a-dependence 
of ^ must differ from the a-dependence of m + by only a multiplicative constant. 
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Three-term 


Requirements for 


Compatible with 


set 


irrediiril)]e T-indenendenre 


a new solution^ 

CAi XXV_j V V k.) \_7 X l_X U X\_7XX • 




Contradiction! must reduce 


No 


il 2 41 


dontradictioTT must reduce 

X^V-fiX UX CAjVJ-XV^ U X WXX . XXXIJ-OU X vX LA. v_y 


No 


11 2 51 


?/ = no, -\- V 


No 


{1,2,6} 


^ = f + 7; 


If a 7^0 


|1,3,4| 


Contradiction: must reduce 


No 


{1,3,5} 


Contradiction: must reduce 


No 


{1,3,6} 


Contradiction: must reduce 


No 


{1,4,5} 


= + + 7; 

_ . 

da; Xi+/3 ~ r*" 


If a ^ 0; 5 ^ 


{1,4,61 


2 ' ' a 

1 X2+/3 




il 5 61 




No 


(2 3 41 


Contraciictioii' rnust roducf^ 

VyXX UX C_XV_X.X V_y UXVy XX • XXX l_ I- i.f \J X v_ v. L l_ I- 


No 


{2,3,5} 




If a 7^ 0; /3 7^ 0; 




Xn — IIP^^ 4- V 






X? + aXo — a 




-f9 


v^uiiti aLiictiuii. iiiUbL icLiuct; 


No 










_dXi 1 jv" = „ 




{2,4,6} 


^ = f + /3 + ^; 


If a ^ 0; /3 ^ 




Xo = ue^'' + z/- 


u ^ 




_dX2 I /^J^ ^ ^ 




|2,5,6| 


= aX^ + 5 


No 


{3,4,5} 


Contradiction: must reduce 


No 


{3,4,6} 


M = + ^ 


If a 7^0 


{3,5,6} 


Contradiction: must reduce 


No 


{4,5,6} 


= aXi + /3 


No 



Table 4. Assuming that a three-term set is irreducibly independent of x results in 
restrictions on the possible solutions. The first column lists each possible three-term 
set. The middle column shows the consequences necessary for that set to be irreducibly 
independent of a;. If this requirement is compatible with solutions not included as one 
of the special cases, restrictions on such a solution are shown in the third column. 
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in which case X2 + aXi = (3 for some constants a and (3. By shifting the zero of a, we 
absorb X2 into Xi by shifting the zero of a. So this particular case can be reduced to 
the special case of X2 = and does not yield any new solutions. 

By examining all possible four-term sets, we find that no new solutions are admitted 
by any combinations that are not included as one of our special cases. Applying the 
same method to five-term sets again yields no new results. Thus, we have examined all 
possible combinations of one-, two-, three-, four-, five-, and six-term sets and discovered 
that they allow no solutions other than those covered by one of the special cases. Since 
the special cases include all known conventional potentials, we therefore conclude that 
this method finds all of such known potentials, and it precludes other solutions. We 
have thus proven that the set of known /i-independent solutions is complete. 



4. Superpotentials that contain explicit h dependence 

Thus far we have found all known additive shape-invariant superpotentials that do not 
depend explicitly on and have proven that no more can exist. We now show that 
our formalism can be generalized to include "extended" superpotentials that contain % 

explicitly. In this case, we expand the superpotential W in powers of %: 

00 

W{x, a,h) = Y^ h'^Wnix, a) . (18) 

n=0 

We wish to substitute (HID in (jH]). From ([H]), 

dW _^ dWnjx.ao) 
dx ''^='^" dx 

and 

00 00 

{x,ao,h)=Y.Y.^'^'WkWi. 

1=0 k=0 

Since ai = ao + h, W {x, ai,h) = W (x, + %,%). Expanding in powers of ^, 



W{x,a^,h)= EEtT f),k 

m=Ok=0 - 



a=ao ' 



So 



Similarly, 



00 n s tn f)^~^ 

(X, a„ = E E E i^^^s-k) 



dx ±^at-1. k\ da^dx ' 



After significant algebraic manipulation and requiring that the result must hold for 
any value of h, we find that the following equation must be true separately for each 
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positive integer value of n: 

dWn-l 



k=0 



dx 




+ £ W^.a^. - I :^ 3i I = o- (19) 



For n = 1, we obtain 

- 1; K 

yielding '^ q%^\'q^ = §^ + g) for k > 1. We have shown that all conventional 
superpotentials W = Wq are solutions of this equation. The extended cases (PE]) are 
solutions to f lT9|) as well. In addition, new potentials can be generated by applying f|T9l) 
for all n > 1, as we show in the next section. 

5. Generating a new h dependent superpotential 

We now use the method outlined in the previous section to generate a new potential. 
We begin by choosing a conventional shape invariant solution to satisfy (l20l) : Wq = 
— atanhx — 5secha;, which is a Scarf II superpotential from Table [1] with A = —a. For 
n = 2, the expansion in (fT9|) yields 
dWi d 

and for n = 3, we obtain 

dW2 _ d {2W0W2 + W^) _ 1 d'^WoWi 2 O^Wq _ ^ 
dx da 2 da"^ 3 da^dx 

Without imposing boundary conditions, there are many solutions to f ll9l) for each value 
of n. However, from physical considerations we require that 1) solutions should not 
have singularities worse that 1/x^ to prevent domain splitting or particles being sucked 
into the singularity; and 2) the asymptotic limits of W be the same as those for the 
corresponding Wq, so that supersymmetry remains unbroken. 

With these considerations, these two coupled equations are solved by Wi = 
and W2 = — (^a^B sinhx)^ ' '^^^ next order equations are solved by W3 = and 
Wi = — ^(^a-B'sinhx)^ ' Generalizing these, we get Wq = — atanhx — _Bsechx; and for 
n > 1, 

4_B cosh X 

W2n-i = 0; = " (2a - 25 siuh x^^^ ' 
yielding a sum that converges to 
W{x, a, h) = —a tanh x — i?sech x 

-2Bhcoshxl— ^ , , T-—^ r.^ . X t) ■ (21) 

\2{a — B siim X) — n 2[a — B smhx) + n J 

This is a hitherto undiscovered superpotential that meets the requirements of shape- 
invariance. 
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6. Conclusions 

We have transformed the condition for additive shape invariance into a set of local 
partial differential equations. For conventional cases that do not depend on ft, we have 
shown that the shape invariance condition is equivalent to an Euler equation expressing 
momentum conservation for fluids and an equation of constraint. Solving these equations 
we have generated all known ^-independent shape invariant superpotentials, and we have 
also shown that there are no others. 

For extended cases in which the superpotential depends explicitly on we have 
developed an algorithm that is satisfied by all additive shape invariant superpotentials. 
Finally, we have generated a new shape-invariant superpotential using this algorithm. 
This method thus has the possibility to greatly expand our ability to identify shape- 
invariant superpotentials. 

It may also be possible to extend this method to other forms of shape invariance 
such as multiplicative or cyclic. For these types of shape invariance, the potentials are 
generally not available in terms of known functions, except in very special cases {N — 2 
for cyclic and limiting cases for multiplicative). It remains to be shown whether for 
these, the shape invariance condition can be transformed from a difference-differential 
equation into a set of partial differential equations and be subjected to similar analysis. 
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